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Abstract 

We present an algebraic Bethe ansatz for certain submanifolds of the bilayer vertex 
models proposed by Shiroishi and Wadati as coupled asymmetric six-vertex free-fermion 
models. A peculiar feature of our formulation is the presence of a diagonal monodromy 
matrix element that does not generate unwanted terms. The model contains two free- 
parameters entering into the Bethe ansatz equations as a pure phase factor. 
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The six- vertex model satisfying the free-fermion condition possesses certain special proper- 
ties which are not common to general integrable vertex models of statistical mechanics [1] . This 
model provides a three-parameter family of commuting transfer matrices [2, 3] whose Boltz- 
mann weights satisfy generalized forms of the star-triangle relations [4, 5]. The later property 
is the basis to generate new exactly solvable vertex models by coupling in an appropriate way 
a pair of six- vertex free-fermion models [4, 5]. This remarkable possibility has originally been 
formulated by Shastry [4] in the case of totally symmetric six-vertex free-fermion model. As 
a consequence of that Shastry was able to rederive the i?-matrix of the covering vertex model 
for the one- dimensional Hubbard chain from a rather promissing perspective [4] . 

More recently, Shiroishi and Wadati [6] extended such an approach to include the asymmet- 
ric six-vertex free-fermion model which have resulted in new interacting bilayer vertex systems. 
To our knowledge, however, the Bethe ansatz solution of such potentially interesting integrable 
bilayer systems is still an open problem. The purpose of this letter is to make a first step toward 
solving this problem by formulating the algebraic Bethe ansatz for what we expect to be the 
simplest branches of the bilayer vertex models discovered by Shiroishi and Wadati [6]. These 
submanifolds are interesting both because it does not contain the covering Hubbard model 
and due to an unusual feature concerning the quantum inverse scattering method. It turns out 
that one of the diagonal elements of the monodromy matrix does not produce unwanted terms 
thanks to a remarkable annihilation property of the eigenvectors. The Bethe ansatz equations 
are of rational type with an additional phase-shift that encodes all possible free-parameters 
dependence. 

The building blocks of the bilayer vertex model constructed by Shiroishi and Wadati are 
two trigonometric Felderhof 's free-fermion model [2] whose Lax operators are defined by 



on the site i of a lattice of size N. Here A denotes the auxihary space, A is the spectral 
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parameter and 7^ play the role of color variables that are attached to each k-th site of the lattice. 
The functions u;±(A, 7^, 7,-) = °+(^.7..Ty)±''+(^.7..Ty) and a;±(A, 7,, 7,) = fe-(A.7.,7.)±a-(A,7.,7,) ^^^^^ 

c^ilA, 7j, 7j) and 6±(A, 7^, 7^) are the elementary Boltzmann weights parameterized in terms of 
the spectral parameter A and the color variables 7^ by 

a±(A, 7j, 7j) = cosh(7i - 7^) cos(A) ± sinh(7j + 7^-) sin(A) (2) 

^±(A, 7i, 7j) = cosh(7i + 7j) sin(A) ^ sinh(7j - 7^) cos(A) (3) 

The Lax operator of the bilayer vertex model is given by a suitable coupling of two Felder- 
hof 's free-fermion models and it can be written as follows 

^AiiK, k+i, li, li+i) = -^2(Ai - Aj+i, 7i, 7i+i)L5i](Ai - Aj+i, 7i, 7^+1) 

(Ai + Ai+i, 7i, -7i+i)(jf _4L5[](Ai + Ai+i, 7i, -7i+i)(7j^_4 
+h[L^Ai{K + Aj+i, 7i, -7i+i)(Tf^_4L5i](Aj - A^+i, 7^, 7^+1) 

(Ai - Ai+i, 7i, 7i+i)L5j[](Ai + Ai+i, 7i, -7i+i)o-i,^] (4) 

The solution found by Shiroishi and Wadati [6] for the couplings /i 2 depends on two arbi- 
trary constants denoted by A and B in their work. Here we are going to consider the simplest 
submanifolds as far as the functional behavior of couplings is concerned. They occur when 
one of the constants is null because the other can be rescaled to unity without loosing gener- 
ality. Even in this situation there exists at least two free-parameters since the Lax operator 
(4) depends on both the difference and the sum of the spectral parameters and color variables. 
It is convenient to define the couplings of these submanifolds with the help of the relations 
fl'±(A, 7) = 0^(^70) • ■'■^ terms of these auxihary functions they are given by 

f ^[^±%l^ + ^^%l-^-2] forB=0 

i[^±i^ll^ + ^^i^^-2] forA=0 
I 4 5(+(Ai+i,7i+i) 5(_(Ai+i,7i+i) 

f l[^+(Am,7m)_^;-(Am,7m)j 
/■ _ J 4 .(?+(A,,:,7,) ^-(A,-.7i) 

I 4 c/_(Ai+i,7j+i) t/+(Ai+i,7j+i) 
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and 



The fact that the L-operator (4) satisfies a colored version of the Yang-Baxter equation [7] 
imphes that it is possible to define a family of row-to-row commuting transfer matrices T{X) 
as the trace of the monodromy matrix T(A) on the auxiliary space of the coupled six vertex 
free-fermion models. The monodromy matrix then reads 

T(A) = Cal{\ Ao, 7, 7)^AL-i{K Ao, 7, 7) ■ ■ ■ ^Ai{K Aq, 7, 7) (7) 

These square lattice vertex models have two free parameters Aq and 7 and their integrability 
is guaranteed by the following intertwining relation 

R{\, ij)T{X) ® T(/i) = T(/x) ® T{i^)R{X, /x) (8) 

where the i?-matrix is R{X, //) = Pi2>Ci2(A, 7, 7) and P12 is the permutation operator. 

We now turn to the algebraic Bethe ansatz solution of these vertex models. The tech- 
nicalities entering the solution of the submanifolds ^4 = and B = are very similar and 
these details will be presented only for the former case. The main results, however, such as 
the eigenvalues of the transfer matrices and corresponding Bethe equations will be given for 
both of them. An essential ingredient in this approach is the existence of a reference state in 
which the monodromy matrix acts triangularly [8, 9]. This helps us to identify the elements of 
the monodromy matrix as particle creation and annihilation fields. Here we denote this state 
by |0) and choose it as the standard ferromagnetic vacuum where all the spins are in the up 
eigenstate of o"^ j. The action of the monodromy matrix in this state yields the relation 

/k(Ar t t t \ 

[u;2{X)f t 
[u;2{X)f t 

V [cos{X)]^J 

where the symbol f denotes non-null values and the functions a;i(A), a;2(A) and a;3(A) are given 

by 

uJi{X) = a+(A - Ao, 7, 7) + a+(A + Aq, 7, -7)[/i«+(A + Aq, 7, -7) + 2/2a+(A - Aq, 7, 7)] (10) 



T(A) |0) 



|0) (9) 
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uj2{X) = a+(A- Ao,7,7)^'-(A- Ao,7,7) + /ia+(A + Ao,7, -7)^'-(A + Ao,7, -7) 

/2[a+(A + Ao, 7, -l)b^{X - Ao, 7, 7) + a+(A - Aq, 7, l)b-{X + Aq, 7, -7)] 

(11) 

a;3(A) = bl (A - Ao, 7, 7) + (A + Ao, 7, -7) [fib- (A + Ao, 7, -7) + 2/26- (A - Ao, 7, 7)] (12) 

The next step is to seek for a suitable representation of the monodromy matrix. Previous 
experience in deahng with the quantum inverse scattering method for genuine four dimensional 
representation [10] suggests us to take the following form 

/5(A) B{X) F{X) \ 
r(A) = C{X) i(A) B*{X) (13) 
U(A) C*{X) D{X)),,, 

where as in ref.[10] -B(A), B*{X), C(A) and C*(A) are two dimensional vectors, the field A{X) 
is a 2 X 2 matrix and the remaining operators are scalars. 

Comparison between the above representation and the triangular property (9) reveals that 
B{X), B*{X) and F{X) play the role of creation operators with respect to our choice of ref- 
erence state. This by no means is an indication that all of them are needed to generate the 
eigenstates of the transfer matrix. To make further progress we have to investigate the com- 
mutation relations between the monodromy operators from the Yang-Baxter algebra (7). In 

order to present the results as general as possible we find convenient to represent the i?-matrix 

4 

by i?(A, /x) = ^ i?^^(A, /x)e(jc ® where Cab are the Weyl matrices. We first look at the 

a,b,c,d 

commutation rules between the creation operator B(X) with itself and with the diagonal fields 

and the results are 

B{X)^B{f,) = ^^^[B{^)^B{X)].f{X,^i) (14) 

^^^^^^^^ ^ mt^f^''^^^^^ " m^f^^^^^''^ ^'^^ 
^^^^ ® ^^^^ = M^)^^^"^®^^^^^-'^^'^^ " muiff^^^®^^"^ ^^^^ 
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r{X,IJ,) 



The auxiliary matrix r(A, /i) is given by 

/I 0\ 

a(A,/i) b{\, fi) 

b{X,fi) a{X,fi) 

Vo ly 

whose non-trivial matrix elements a{X, /i) and b{X, /i) are 

cos(2A) — cos(2/i) 



(18) 



a(A,//) = 



cosh (27) 



&(A,//) = 



(19) 



+ cos (2 A) - cos(2/x) + cos (2 A) - cos(2/x) 

Prom these results we first recognize that the matrix r(A, n) is related to the one of the 
isotropic six-vertex model via reparametrization of the weights in terms of a new variable 
A = cos(2A). Next wc observe from the commutation rules (14-16) that two creation fields 
B{X) can not generate other creations fields and also that all the unwanted terms coming from 
the diagonal fields -B(A) and ^(A) remain in the space of the states of the operator -B(A). 
It is notable that such commutation relations are in the same form of that appearing in the 
nested Bethe ansatz solutions of multistate generahzations of the six- vertex model [11, 12]. A 
relevant difference, however, is the presence of the diagonal field -D(A) which requires further 
elaboration. As far as the one-particle state is concerned, the commutation relation (16) 
shows that this operator generates unwanted terms that are trivially canceled out due to the 

— * 

annihilation property C*(A) |0) = 0. To show that this is indeed a general situation we need 
the help of an extra commutation rule, namely 

c*ix) B{^.) = ® c*W - ® (20) 

Now by iterating this relation on the tensor product of several B{Xi) we are able to derive 
the following remarkable annihilation property 



C*(A) (8) B{Xi) ® . . . (8) B{Xn) |0) = 



(21) 



which together with the commutation rule (17) implies that the field D{X) does not produce 

— * 

unwanted terms in the space of states generated by the creation fields B{X). 



These observations strongly suggest us to take as the eigenvectors of the transfer matrix 
T(A) the following hnear combination 

|$(Ai, . . . , A„)) = S(Ai) (8) B{X2) <g) . . . B{X^).f |0) (22) 

where the vector ^ depends on the parameter Aj and its components will be denoted by ^"i """ 
where each = 1,2. 

At this point we have gathered the basic tools to determine the eigenvalues A(A, {A^}) of 
the transfer matrix defined by 

[B{X) + A„„(A) + D{X)] |$(Ai, . . . , AJ) = A(A, {A,}) |$(Ai, . . . , AJ) (23) 

a=l 

In order to do that we use the commutation rules (15-17) and the property (21) to commute 
the diagonal operator through all the B{Xi) fields until they hit the reference state. Their values 
acting on the reference state can be read out from Eq.(9). Omitting details of the calculations 
our final result for the eigenvalues of the transfer matrix is 

A(A, {A,}) = [.,{xr n §rT^,+^-^^^r n llrYi+t'^^^)]" n llrTi^'''^^' 

(24) 

provided that the rapidities satisfy the Bethe equations 

=(-1)"-^AW(A = A„{A,}), z = l,...,n (25) 

The function A'^^)(A, {A;}) refers to the eigenvalue of an auxiliary problem related to the 
condition that the vector be an eigenvector of the transfer matrix whose Boltzmann weights 
are the i?-matrices (18), namely 

rl\l{\ Ai)<-^(A, X,) . . .<7r"(A, \n)r"' = A(^)(A, {A,}).^''"-^^ (26) 

This problem is solved by another Bethe ansatz which gives the well known results [8, 9] 
for the inhomogeneous rational six-vertex model, 

m -| n m -I 

A«(MM. w) ^ n^+nKA.A.,n^ 127, 
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where the auxihary variables {/ij} are subject to the following set of constraints 



nK«.A.) = -nS^, , = 1. 



, m 



(28) 



i=i fc=i 

Now we first insert the nested eigenvalues (27) into Eq.(24). By taking into account the 
explicit expression of the i?-matrix elements after carrying out cumbersome simplifications we 
find that 



A(A, {A,}, {M) = MX)]'^l[z{Xn) 



U - cos(2Ai) + cos(2A) 



+MA)]^<ln^(A,7) 



cos(2Aj) — cos(2A) 

+ cos(2A) - cos(2Ai) 



+n-^(A'7)n 
j=i 1=1 



cos(2A) — cos(2A 
?7 + cos(2A) -cos(2/iz)" 



n 



U + cos(2/x,) - cos(2A) 



f-J^ cos(2A) — cos(2/x;) 



cos(2/x/) — cos(2A) 

where the function z{X,'y) and the parameter U are given by 
6+(A,7,0) 

^(A,7) 



(29) 



a_(A,7,0) 
a+(A,7,0) 



for B=0 



and U — < 



for A=0 



cosh (27) 
2 

cosh(27) 



for B=0 
for A=0 



(30) 



6_(A,7,0) 

Proceeding similarly for Eqs. (25,28) we obtain the corresponding nested Bethe ansatz equa- 
tions 

1 N 



U — cos(2Aj) + cos(2Ao 
— cos(2Aj) + cos(2Ao) 



[ .(\ ( 1 ^n-l + cos(2/x;) - cos(2A0 



" cos(2Aj) - cos(2/Xi) 
M U + cos(2//,) - cos(2Aj) 



cos(2Aj) - cos(2^,) 

^ C0S(2/X;) - COs(2/Xfc) - U 

cos(2///) - cos(2//fc) + [/ ' 



(31) 



Prom the result (31) we observe that it is possible to cast the Bethe equations in a simpler 
form if we define the following new rapidities 



cos(2Aj) ~ 1 cos(2Ao) cos(2/ij) _ cos(2Ao) 
U ^ '^2^ U ' U ^^'^ U 



(32) 
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leading us to the expressions 



N 



l=\ K — p'l + \ 




[-^(Ao,7o)] 



N 




n 



(-irn 




k=l 



n 



m 



jll- jlk + 1 



I — 1, ... ,m 



(33) 



In terms of these variables the dependence of the Bethe equations on the free-parameters 
becomes restricted as a fixed phase factor. Next we observe that their rational form resemble 
that of the supersymmetric t-J model [13, 14] with the important difference that in our case 
we have a bigger space of states such that m < N rather than m < N/2. As far as the 
Hilbert space is concerned one would be tempted to compare these equations with that of the 
supersymmetric U model [15] . However, this demands that the continuous parameter governing 
the four dimensional representation of the superalgebra spl{2\l) reaches the value related to 
atypical representations which is not permitted in this model [16, 17]. These remarks are 
strong evidences that the class of models solved here lies in between these two supersymmetric 
systems motivating subsequent investigations. One of them consists in the computation of 
physical properties of the bilayer vertex model such as the ground state behavior dependence 
on the free-parameters. Another interesting point raised here is the possible connection between 
this system and four dimensional atypical spl{2\l) representations. It seems also worthwhile to 
generalize our results to include the manifold discovered by Shiroishi and Wadati for arbitrary 
values of the constants A and B. It is plausible to beheve that this will lead us to Bethe 
equations interpolating between that of the t-J and Hubbard like models. In this situation a 
crucial point is to unveil the corresponding six-vertex hidden symmetry which has eluded us 
so far. 

Acknowledgements 

This work has been partially supported by the Brazilian research Agencies CNPq and Fapcsp. 



8 



References 

[1] C. Fan and F.Y. Wu, Phys.Rev. 179 (1969) 560 
[2] B.U. Felderhof, Physica 66 (1973) 279, 509 

[3] V.V. Bazhanov and Y.G. Stroganov, Theor.Math.Fiz. 62 (1995) 337 
[4] B.S. Shastry, J.Stat.Phys. 30 (1988) 57 

[5] I.G. Korepanov, Mod.Phys.Lett.B 3 (1989) 201, Comm.Math.Phys. 154 (1993) 85 
[6] M. Shiroishi and M. Wadati, J.Phys.Soc.Jpn. 64 (1995) 4598 

[7] J. Murakami, Int.J.Mod.Phys. A (1992) 765 ; T. Deguchi and Y. Akutsu, J.Phys.Soc.Jpn. 
60 (1991) 4051 

[8] L.A. Takhtajan and L.D. Faddeev, Russ.Math.Surv. 34 (1979) 11 

[9] V.E. Korepin, G. Izergin and N.M. Bogoliubov, Quantum Inverse Scattering Method, Cor- 
relation Functions and Algebraic Bethe ansatz, Cambridge University Press, Cambridge, 
1992 

[10] P.B. Ramos and M.J. Martins J.Phys.A-.Math.Gen. 30 (1997) L195 

[11] P.P. Kulish and N.Yu. Reshetikhin, J.Phys.A: Math.Gen. 16 (1983) L591 

[12] O. Babelon, H.J.de Vega and CM. Viallet Nucl.Phys.B 200 (1982) 266 

[13] B. Sutherland, Phys.Rev.B 12 (1995) 3795; P. Schlottmann, Int.J.Mod.Phys. B 11 (1997) 
355 

[14] F.H.L. Essler and V.E. Korepin, Phys.Rev.B 46 (1992) 9147, A. Foerster and M. Karowski, 
Phys.Rev.B 46 (1992) 9234 

[15] A.J. Bracken, M.D. Gould, J.R. Links and Y.-Z. Zhang, Phys.Rev.Lett. 74 (1995) 2768; 
R.Z. Bariev, A. Kliimper and J. Zittartz, Europhys.Lett. 32 (1995) 85 

9 



[16] Z. Maassarani, J.Phys.A:Math.Gen. 28 (1995) 1305 

[17] K.E. Hibberd, M.D. Gould and J.R. Links, Phys.Rev.B 54 (1996) 8430; P.B. Ramos and 
M.J. Martins, Nucl.Phys.B 474 (1996) 678; M.P. Pfannmiiller and H. Prahm, Nucl.Phys.B 
479 (1996) 575 



10 



